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The boson peak in metallic glasses is modeled in terms of local structural shear rearrangements. 
Using Eshelby's solution of the corresponding elasticity theory problem (J. D. Eshelby, Proc. Roy. 
Soc. A241, 376 (1957)), one can calculate the saddle point energy of such a structural rearrange- 
ment. The neighbourhood of the saddle point gives rise to soft resonant vibrational modes. One 
can calculate their density, their kinetic energy, their fourth order potential term and their coupling 
to longitudinal and transverse sound waves. 
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I. INTRODUCTION 

As yet there is no generally accepted explanation of 
the boson peak in the neutron or Raman scattering in- 
tensities of glasses. The boson peak is a broad peak at an 
energy transfer of a few meV, where simple crystals show 
only sound waves. Glasses seem to have a sizable excess 
of vibrations at this boson peak. At present, there is no 
agreement which forces drive these extra vibrations into 
the low-frequency region, though several possible expla- 
nations have been proposedii£i2ii££i££. Another contro- 
versial question^iiS is whether the interaction with these 
vibrations is the physical reason for the Ioffc-Regel limit, 
the reduction of the mean free path down to the wave- 
length of sound waves in the THz range. 

The present paper proposes a detailed atomic model 
for the boson peak modes in glasses consisting of close- 
packed atoms or spherical molecules: The soft modes are 
ascribed to small regions with a pronounced shear misfit 
with respect to the surrounding matrix. In the following 
section II we describe the picture and derive the prop- 
erties of these modes. Section III treats the connection 
to the extended soft potential modeli and compares the 
predicted sound wave scattering and softening to exper- 
iment. Section IV discusses and summarizes the results. 



position after a distance of 2r/v3 (r sphere radius), the 
situation at the right side of Fig. 1. The six atoms there 
form a pair of edge-sharing tetrahedra. Since the dis- 
tance between close-packed planes is 2ry/2/3, the shear 
angle is 1/V2 in radian, an angle of 40.5 degrees. 

The physical problem of a small piece of matter able 
to transform to a sheared shape within an elastic matrix 
has been treated fifty years ago by several authors, no- 
tably by J. D. Eshelb y 11 ! 12 . Here, we translate Eshelby's 
result into the usual convention, in which the shear angle 
e and the shear stress a are related by a = Ge (G infi- 
nite frequency shear modulus). Let Vi be the volume of 
the spherical inclusion and e% the shear angle difference 
between its two stable configurations (in the example of 
Fig. 1, a — l/v2 and Vi — 6v, where v is the atomic 
volume). For strictly harmonic potentials, the symmetric 
case of two equally strained configurations has the energy 



E a = lGv. t ej 

The coefficient 7 is given by 
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where op is Poisson's ratio. Since Poisson's ratio lies 



II. THE GLIDING-TRIANGLE MECHANISM 
A. Shear strain defects 

The central concept of this model is a structural re- 
arrangement of a limited region in the sample which 
changes its shape to a sheared one. In three-dimensional 
close packing, the lowest-energy structural rearrange- 
ment can be viewed as a gliding of a triangle of three 
close-packed atoms over an underlying close-packed plane 
(Fig. 1). On the left side of Fig. 1, the center of the up- 
per triangle lies exactly over the center of a lower triangle 
within the close-packed plane. The central points of the 
six atoms form an octahedron. Gliding the triangle into 
the vertical direction in Fig. 1, one reaches again a stable 




FIG. 1: The gliding-triangle rearrangement of six closely 
packed spherical atoms. 
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between 0.1 and 0.44 for the known glasses, 7 lies between 
0.48 and 0.57, close to 1/2. 

Eshelby's solution divides the energy into two almost 
equal parts, one located in the inclusion and one outside. 
Their ratio is 7/(1 — 7). In the symmetric case, the in- 
clusion would have to distort by e,/2 to fit exactly into 
the unstrained hole in the surrounding matrix. From a 
structural point of view, this is the saddle point between 
the two stable configurations. 

In strictly harmonic approximation, Eq. JT]) holds and 
the structural saddle point indeed is a saddle point of the 
energy landscape. However, the harmonic approximation 
does not account for the flattening at the saddle point, at 
least as regards the distortion energy E c of the inclusion. 
For the inclusion, a better description is by a cosine law 
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shown in Fig. 2, where e c is the shear angle of the central 
inclusion, counted from the saddle point of the energy of 
the inclusion. This cosine law has the correct second 
derivative of 6Gv at e c = ±1/2-^/2, the two structural 
energy minima of Fig. 1. 

The embedding energy outside the inclusion depends 
on the difference between e c and the zero point e m of the 
shear of the matrix. In terms of e c , the shear state of the 
hole after removing the inclusion and relaxing the matrix 
is e c = e n . Let us assume 7 = 1/2. Then the second 
derivative of this outside embedding term must also be 
6Gv. The total elastic energy E e i is given by the sum of 
inside and outside contributions 
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Note that the symmetric case at e m — (the red lines 
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FIG. 2: Potential energies for inclusion and matrix as a func- 
tion of the inclusion shear e c around the structural saddle 
point, where the matrix equilibrium shear e m lies in the mid- 
dle between the octahedron and the two tetrahedra of Fig. 1. 
The continuous black line is the inclusion energy, the dashed 
red line and the dashed blue line show the matrix energy for 
e m — and 0.05, respectively. The red dot-dashed and blue 
dotted lines show the corresponding total energies. 



in Fig. 2) has its minimum at e c = 0, with all the en- 
ergy inside the inclusion and a vanishing restoring force 
constant for a displacement of e c . Thus one does not 
get a saddle point of the energy landscape, but one gets 
a shear mode with zero restoring force constant at the 
structural saddle point. In this case, there is no elastic 
energy outside the inclusion and the creation energy E s 
for the soft mode configuration is 



E„ = 



3Gv 
2?r 2 ' 
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In order to know the probability of finding such a soft 
gliding-triangle mode in a closed packed glass, one needs 
to know the ratio Gv/ksTg, where T g is the glass tran- 
sition temperature, at which the glass falls out of the 
thermal equilibrium. For metallic glasses, this informa- 
tion can be extracted from a recent data collection^. 
Fig. 3 plots the product Gv for 30 metallic glasses at the 
glass transition as a function of T g . There is a marked 
scatter, but on average Gv/ksTg — 17.6 (the line in Fig. 
3). With eq. ©, this implies a ratio E s /k B T g of 2.67, 
i.e. a Boltzmann factor of 0.069. Thus one has to reckon 
with a sizeable number of nearly unstable resonant shear 
modes in a metallic glass. 



B. Properties of soft shear modes 

The equilibrium value e c o of the inner coordinate e c 
for a given e m o is determined setting the first derivative 
of Eq. (0J with respect to e c equal to zero 



e m o — e c o 
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The second derivative of Eq. (|4]) at e c rj gives the restor- 
ing force constant 
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FIG. 3: Energy Gv as a function of the glass transition tem- 
perature T g in metallic glasses^. The line depicts the average 
Gv/ksTg = 17.6. 



3 



which vanishes for e m o = e c o = 0. 

To calculate the kinetic energy of the mode, consider 
a sphere with a volume of 6u, corresponding to the six 
atoms in Fig. 1. Its radius tq is given by 



4-7rrj] 



6v. 



(8) 



Inside the sphere, the motion is a pure shear with am- 
plitude e = e c — e c o, say u x — ey/2, u y — ex/2 and 
u z = 0. The motion of the outside has an amplitude de- 
caying with the square of the distance from the center, 
which joins the inner shear amplitude continuously at the 
sphere boundary. The kinetic energy inside the sphere is 
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where p is the density. The kinetic energy outside the 
sphere, in which the r 4 of Eq. (fT0|) is replaced by Tq / r 2 
and r is taken from tq to infinity, is a factor of five higher 
than the one inside. Replacing the density p — M/v, 
where M is the atomic mass, one finally arrives at 
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For many purposes, it is useful to define a soft mode co- 
ordinate A with the normal coordinate property Ekm = 
A 2 j 2. For the soft shear modes 
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Consider the shear mode at the structural energy min- 
ima, say the octahedron of Fig. 1 at e m = —l/2\/2. 
Then, D e = 12Gv, and the mode frequency attains its 
maximum uj mnT with 
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It is interesting to compare this frequency with the 
Debye frequency lud, determined essentially by the trans- 
verse sound modes. The transverse sound velocity v t is 
given by 



2 G Gv 
Vt = ~p = M ■ 
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The Debye sound velocity vd is given by the average 
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where vi is the longitudinal sound velocity, on the average 
for the metallic glasses 13 a factor of 2.26 higher than the 
transverse one. For this case, vp = 1.13v t . The Debye 
wave vector is ko = (67r 2 /v) 1 / 3 , so on average 



2.72 uj n 



(15) 



This shows that the gliding-triangle mode considered 
here has a maximum frequency of less than half the De- 
bye frequency in the absence of any structural strain. 



C. Density of gliding-triangle modes 

In order to estimate the density of soft shear modes, 
one first needs the number of octahedra per atom in 
the disordered structure. In the crystalline close-packed 
structures, both fee and hep, one has one octahedron per 
atom. Though this number might easily be smaller in 
the disordered case, we will take it there also as one. 

For each octahedron, there are twelve possibilities to 
distort to a pair of edge-sharing tetrahedra. There are 
four choices of the two triangles which are to glide against 
each other and for each pair there are three glide direc- 
tions, 120° rotated to each other. Consequently, there 
are two stable minima per atom, the octahedron and the 
edge-sharing pair of tetrahedra, and twelve saddle points 
between them. As shown in the Appendix, the twelve 
saddle points lie pairwise rather close to each other, with 
an angle of only 38.9 degrees between the two members 
of each pair. 

At the glass temperature T g , the saddle points are al- 
ready sufficiently high in energy to be, in first order, neg- 
ligible for the partition function, which can be approx- 
imated by a harmonic potential around the octahedron 
energy minimum at e oct — —1/2^/2. Close to the min- 
imum, Eq. © gives e c0 - e oct = (e m0 - e oct )/2. At 
the energy minimum, Eq. ^ takes the simple quadratic 
form 
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with the partition function 



Z = 



27rk B T g 
3GV 



(16) 



(17) 



at the glass temperature T g . 

The probability p s — p(e m — 0) to find e m close to the 
structural saddle point is 



3exp(-E s /k B T g ) / 54Gv 



■ exp 
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The factor 3 stems from the fact that there are three 
saddle point pairs per minimum, as explained in detail 
in the Appendix. 

It is useful to introduce the coordinate x 



x = \[2'Ke C Q 



(19) 



which is zero at the structural saddle point and ±7r/2 at 
the two energy minima. From equs. (JT]) and (TT2]) . one 
obtains for the frequency of the shear mode 



W = 0J max SIM. 



(20) 



As one moves away in e m from the structural saddle 
point, the probability density changes only slowly, but 
x and the restoring force constant of the mode change 
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rapidly (see Fig. 2). 
related to the one of e„ 
lated from eq. (|6|) 



The probability density of x is 
by the derivative de m /dx calcu- 



dx 
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The vibrational density of states of the soft shear 
modes is given by 



p s de m0 dx p s uj 2 
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where the prefactor 1/3 results from the usual convention 
of normalizing the three modes per atom to 1. Since we 
are only interested in the density of states well below 
u> m , we neglect the square root term. Then g 8 (uj)/uj 2 is 
simply constant 



Ps 



(23) 
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This constant is to be compared with the constant 3/q 
obtained by dividing the Debye density of states by to 2 . 
For the average Gv/ksTg = 17.6 (see Fig. 2), one finds 
with Eqs. (fl~5|) and (| 18[) a ratio of 0.6 between the density 
of the soft shear modes and the Debye one. Thus one 
has a rather large number of soft gliding triangle modes 
at low frequency, comparable with the number of sound 
waves. 



D. The shape of the boson peak 

The soft shear modes are not exact eigenmodes of the 
harmonic vibrations of the glass. They interact with the 
other vibrations and in particular with the sound waves. 
The low frequency shear modes can be understood as the 
cores of quasi-localized or resonant vibrations. They are 
bilinear ly coupled to the other modes. This system of 
localized soft vibrations coupled bilinearly to extended 
modes (sound waves) is the basis of the soft potential 
modelSi^^. 

Through their interaction with the sound waves, there 
is also an interaction between the soft shear modes. The 
bilinear interaction of the soft modes with the much 
larger number of higher frequency modes causes a down- 
shift of the soft modes, some of which even become har- 
monically unstable. Stabilization via anharmonicity re- 
sults in a soft mode spectrum linear in u> below a fre- 
quency lo c given by the interaction strength, while at low 
frequency one has a DOS oc w 4 . The crossover between 
the oj 4 and Li regions then leads to the boson peaki^l 
The total number of soft modes is not changed by these 
frequency shifts. It has been shown that this interaction 
mechanism can explain the nearly universal density of 
observed two level systems in glasses^ 

For the purpose of the present paper, it suffices to know 
that the exact shape of the boson peak requires a more 



detailed theoretical treatment. Here, we restrict our- 
selves to the conclusions which are easily accessible from 
our postulate for the eigenvector of the excess modes. 
As will be seen in the next subsection, these conclusions 
include a calculation of the damping and softening of 
the sound waves from the measured spectrum of a given 
metallic glass. 



E. Coupling, damping and softening 

The coupling of a gliding-triangle soft shear mode to 
an external shear strain is given by the derivative of the 
elastic energy of eq. (UJ) with respect to both the external 
strain em and the internal coordinate e r 
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if the external shear strain happens to lie exactly in the 
direction of e c . In the general case, the directional aver- 
age over the five possible shear strain orientations must 
naturally be taken into account. 

It is usual to express the coupling in terms of a product 
A t Ae m of the shear strain and the soft mode coordinate 
A, defined via eq. (fTT|) . From A t , one can define a fre- 
quency u> t via 



U3, 



_Af_ 
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Sorting all the factors out, one finally gets 



u>t = 
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about a factor of 3.5 lower than the Debye frequency. 

With u>t, the damping T of the transverse sound waves 
(the full width at half maximum in a Brillouin experi- 
ment) is given by 



r t = 



37T 



Wt 2 5 s (w), 



(27) 



so the treatment predicts a damping of the transverse 
sound which increases with the frequency or wavevector 
squared. The Ioffe-Regel condition r = oj/it is reached 
at the frequency LOjRt 



9TT 2 auj 2 



(28) 



Here a is the ratio between Debye density of states and 
soft shear wave density of states, which in the average 
metallic glass should be close to 1. For the average case 
with Gv/ksTg = 17.6, ujjRt lies at 0.23 u>d- 

The Ioffe-Regel limit for the longitudinal sound waves 
lies markedly higher, because the gliding-triangle modes 
couple only to the shear (or essentially so; we come back 
to this point in the discussion). Since the longitudinal 
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elastic constant is given by C\\ — B + AG/ 3, where B 
is the bulk modulus, a pure shear coupling implies = 
4Aj/3 for the longitudinal coupling constant. With again 
ujf = Af/Mvf, ujf will be a factor of 4v^/3vf smaller than 
ujt, with vi/vt ~ 2.26 a factor of 0.26 smaller—, which 
makes lujm a factor of nearly 4 higher than lojRt- 

The influence of the soft shear modes on the shear 
modulus is given by 

G 



(29) 



where AG is the difference between G above and below 
the boson peak. 



III. COMPARISON TO EXPERIMENT 

A. Vit-4, a heavily studied metallic glass 

The best candidate for a check of these ideas 
is the metallic glass Vit-4, Zr4 6 . 8 Ti 8 .2Cu7.5NiioBe27.5, 
a glass with an exceptionally low critical cooling 
ratoi^, for which every relevant quantity has been 
measure d 15 ' 16 i 17 ' 18 , though the mechanical^, density^ 5 - 
and neutro n 16 ' 17 data were all taken for slightly different 



compositions. At the glass temperature T g = 615 K, the 
shear modulus is 35 GPa, the bulk modulus 112 GPa^ 
and the atomic volume is 0.017 nm 3 , so fnoo — 27. 5 meV. 
Fig. 4 shows the boson peak at room temperature deter- 
mined from inelastic neutron scattering dat a 16 ' 17 in units 
of the Debye density of states (at room temperature^, 
hu> D = 28.6 meV). 

If one takes the excess over the Debye density of states 
as gliding-triangle modes, one can calculate the corre- 
sponding damping of the sound waves. One finds that 
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FIG. 4: The boson peak of the model glass Vit-4 in compar- 
ison to crystal and liquid, determined from inelastic neutron 
scatteringi^ 1 ^ in units of the glass Debye- level 3/Ep. The 
arrow marks the maximum frequency Umax up to which one 
expects gliding triangle modes. 



the Ioffe-Regel limit is never reached, even for the trans- 
verse sound waves. At the boson peak, the damping is 
about half the Ioffe-Regel limit. Since the longitudinal 
waves are even less damped, one expects rather sharp 
longitudinal Brillouin peaks in an inelastic x-ray exper- 
iment. This was in fact measured for another metallic 
glass 10 , NirsijZi'ciT. In this case, one finds a Q 2 dependence 
above uj max ^ which probably has nothing to do with the 
boson peak. 

Though the boson peak is small, it is crucial for the 
mechanical properties and their temperature dependence 
in the undercooled liquid. One calculates AG/G = 0.08 
from eq. (T2"9"l) with the excess modes from Fig. 4. 
This fits in with the measured difference 18 of 26 % of 
the shear moduli of glass and crystal. The shear mod- 
ulus Griineisen parameter Tq = dhiG/dhiv is found 
to be 6.5 from the temperature dependence of and 
the density^ in the glass phase. Consequently, the 
density ratio 1.027 between crystal and glass explains 
two thirds of the shear modulus difference, the gliding- 
triangle modes explain the rest. 

In the undercooled liquid, the number of gliding- 
triangle modes increases exponentially according to eq. 
(fl~8|) (replacing T g by T). They do not only increase 
their number, but also shift to lower frequency, as seen 
in Fig. 4. However, if one makes a rough estimate by a 
linear extrapolation of the liquid data in Fig. 4 to the 
frequency zero, the increase of the area under the curve 
follows in fact eq. (HJ). Taking AG/G = 0.08 as the ef- 
fect of the gliding-triangle modes frozen in at T g , eq. (fTS)) 
predicts an additional reduction of G of 0.023 GPa/K in 
the melt. This is almost exactly found in experiment 18 
(see Fig. 5). The reduction in the liquid is a factor of 
2.4 larger than the Griineisen expectation on the basis 
of the thermal expansion 15 . The difference can be frozen 
into the glassy state by quenching. It is obviously due 
to configurational changes of the melt. The finding is of 
great importance for the understanding of liquids. Our 
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FIG. 5: Shear modulus measurement^ in glassy and liquid 
Vit-4. The dotted line corresponds to a Griineisen parameter 
Fq = 6.5 in the glass. 



6 



model gives a quantitative physical explanation for it: it 
is due to shear instabilities, stabilized by the surrounding 
matrix and requiring a formation energy of about three 
times ksT gi which reduce the shear modulus. 

The experiment^ reveals also an additional decrease 
of the bulk modulus with increasing temperature, though 
not quite as pronounced as the one of the shear modu- 
lus. If this is also attributed to our shear defects, they 
must couple to the compression. From the datai^, the 
coupling constant ratio A; / A t of these defects should be 
1.69 rather than the value \/4/3 for a pure shear defect, 
close to the value 1.6 found for the tunneling states in 
glasses^. We will come back to this point in the next 
section. 



B. Soft potential model 

The gliding-triangle model makes some rather detailed 
predictions for the low temperature glass anomalies 20 . 
The easiest way to do this part of the comparison to 
experiment is to look at the relation between the picture 
given here and the soft potential model, the standard 
model in the field for which an extensive comparison with 
experimental data already exist s 21 ' 22 ! 23 . 

The soft potential mode l 21 ' 22 i 23 connects the boson 
peak modes to the tunneling states 2 ^ and to low-barrier 
relaxation. It postulates a constant density of first and 
second order terms D\ and D 2 around the value zero for 
the potential 



Esoft = W(D lXs + D 2 x 2 s + x 4 s ). 



(30) 



The coordinate x s relates again to the coordinate A of the 
soft mode with Eki n = A 2 / 2. Let us denote the fourth or- 
der term in this coordinate by V4A 4 . Then the zero-point 
energy of the purely quartic potential is by definition W, 
given by the equality of kinetic and potential energy at 
A 



W 



K 2 ,4 



(31) 



This relation defines both W and Aq . The kinetic energy 
results from the uncertainty principle at a confinement 
within ±Aq. The energy W separates tunneling states 
below W from vibrational states above. The coordinate 

X s = A/Aq. 

The gliding-triangle modes of the present paper are 
also distributed around a purely quartic potential (the 
continuous red line in Fig. 2). Its fourth order term 
2iv l Gve\ 1 obtained by the fourth derivation of eq. Jl]), 
leads in combination with the definition of A via eq. (jTTJ) 
to 



W 



52/3^10/9 
2 5/9 3 20/9 



(Gvy/Htf/Mv 2 / 3 ) 2 / 3 . (32) 



Eq. (|52")l supplies an average value W/Hb of 4.4 K 
for the thirty metallic glasses of the data collection of 



Johnson and Samwer-1 3 -, a value well within the range of 
those adapted with the soft potential model to a number 
of nonmetallic glasses 22 . It combines a very high energy, 
Gv, coming from the fourth order potential term, with 
the very low kinetic confinement energy h 2 /Mv 2 / 3 . In 
this context, it is interesting to note that the fourth order 
term of the quartic potential in Fig. 2 does not result 
from the short range atomic repulsion, but rather from 
the decrease of the curvature of the inclusion energy with 
increasing e c , an anharmonic property of the saddle point 
of the inclusion energy. 

Let us next turn to the tunneling states. Of course, the 
potential of eq. is in no case a double-well potential 
with a low barrier able to form a tunneling state, because 
the negative force constant of the inclusion can at most 
cancel the positive force constant of the surrounding ma- 
trix. In terms of the soft-potential model, the treatment 
in Section II corresponds to the case where the density 
p{D%) is not constant around zero, but rather a 5-function 
at zero (this explains why one gets an w 2 -rise in the ex- 
cess density of states rather than the w 4 -rise of the soft 
potential model). But it was already pointed out in sec- 
tion II that the interaction between the quasilocalized 
modes is responsible for the shape of the boson peak and 
the appearance of tunneling state o 27 ' 28 . However, such a 
detailed modelling is beyond the scope of the present pa- 
per. Here we make no statement on the number of these 
tunneling states. But we can make a statement on their 
coupling to the sound waves. 

For W — 4.4 K, a tunneling state with a splitting 
around 1 K requires D% = —6. The two minima of the 
double-well potential have a distance \J — 2D 2 in the co- 
ordinate x s of eq. (|30p from each other. Using the defini- 
tions of equs. (|3"Tj) and (fTTj) , one gets the corresponding 
distance in e c . Multiplying this distance with the cou- 
pling factor 6Gv of eq. (j4]) and again averaging over the 
five possible shears, one gets the coupling constant 74 of 
the tunneling state to a transverse sound wave 
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(33) 



For the data collection on metallic glasses 13 this yields 
an average value of 0.49 eV, not too far away from the 
average value 0.39 of direct measurements of jt of the so- 
called low-temperature "tunneling plateau" in the sound 
absorption of 18 different inorganic glasses^ 9 .. In this 
field, scientists have always been wondering what kind 
of mode coordinate would be able to give such a strong 
coupling. Here we have for the first time a detailed an- 
swer to this question, at least for an important subgroup 
of glasses. 

The coupling to longitudinal sound waves is described 
by the analogous coupling constant 7/. From the con- 
sideration at the end of section II on a pure shear de- 
fect like the one proposed here, one would expect a ratio 
li/lt = \/4/3 = 1.155. This does not agree with the 
experimental finding 19 7^/74 = 1.6. As pointed out in 
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the previous subsection, the temperature dependence of 
the elastic moduli in the Vit-4 melti£ suggests a factor of 
1.69, so metallic glasses in this respect do not differ from 
the general case. At present, we see no obvious reason 
for this additional coupling to the compression. 



IV. DISCUSSION AND SUMMARY 

As seen in section III, the Ioffe-Regel limit for longi- 
tudinal sound waves for metallic glasses should be by a 
factor of two to four higher than the one for transverse 
sound waves. In fact, an inelastic x-ray Brillouin experi- 
ment in a metallic glass reveals a relatively low longitu- 
dinal sound wave damping^. 

However, the metallic glass result^ is an exception 
rather than the rule^. In the nonmetallic glasses, one 
usually does find the Ioffe-Regel limit close to the boson 
peak, maybe slightly above the boson peak. How is this 
possible? 

The question is related to three ratios: The ratio of 
the boson peak height to the Debye level, the ratio of the 
coupling constants A 2 /A 2 , which for a pure shear defect 
is 4/3, and the ratio vf/v^. As it turns out^, this last 
ratio is about 5.1 in the metallic glasses, unusually high 
(compare silica, where this ratio is 2.3). This and the 
weakness of the boson peak in many metallic glasses are 
the two reasons for the weak longitudinal damping. A 
simulation^ in a soft sphere glass, where the ratio vf/v% 
was even higher than for real metallic glasses, showed a 
very weak damping of the longitudinal waves. 

There seems to be some as yet unexplained coupling 
of the boson peak modes to the compression. As pointed 
out in section III, the ratio A 2 /A 2 seems to be 1.6 2 = 2.56 
rather than 4/3. This feature is missing in our model. If 
one ascribes it to different packing factors for octahedron 
and tetrahedra pair, these should differ by about 18 %. 
This explanation seems not to be very convincing. We 
rather think that the absence of the compression cou- 
pling in our model comes from the neglect of the atomic 
roughness of the shell surrounding the six central atoms, 
which introduces a random coupling to the compression. 

Despite this small uncertainty, the evidence for the va- 
lidity of the gliding triangle model in section III shows 
that it is able to compete with other explanations of the 
boson pea k 2 >M» 5 & Z £ . There is no adaptable parameter; 
one only needs the shear modulus, the atomic mass and 
the atomic volume. The picture provides a microscopic 
basis for the empirical soft-potential model^ i 21 ' 22 i 23 , not 
for all glasses, but for the subgroup of metallic glasses. 

The question is: To which extent is this answer uni- 
versal? How should one generalize it to other glasses, 
glasses with covalent bonds forming a random network, 
molecular glasses with the additional degree of freedom 
of the molecular rotation, to the random chain structure 
of polymers? The answer to the question requires fur- 
ther studies. It might be an intelligent guess, however, 
to assume the following two general features: 



(i) The boson peak and the tunneling modes are always 
of the same nature, as postulated by the soft potential 
model 

(ii) Boson peak modes and tunneling modes are due 
to small regions with a strong elastic shear misfit, lying 
close to a saddle point of the inclusion energy between 
two stable structural minima of the inclusion. 

The validity of these two assumptions would explain 
the universality of the low temperature anomalies in 
glasses, including those at slightly higher temperatures 
where the boson peak modes dominate. 

V. APPENDIX: THE 12 SADDLE POINTS IN 
CUBIC NOTATION 

In an fee close packed crystal, the six atoms of the 
octahedron have the coordinates (±a/2, 0, 0), (0, ±a/2, 0) 
and (0, 0, ±a/2) in units of the lattice constant a, related 
to the atomic volume by a 3 = 4i>. 

Let us first consider the gliding-triangle motion of the 
two triangles in the (lll)-plane along < 211 > in Fig. 
6. The upper triangle ABC moves opposite to the lower 
triangle DEF. This gliding motion involves not only a 
shear, but also a rotation. To get the rotation-free shear, 
let us define the two unit vectors a\ —< 211 > /V6 and 
a 2 =< 111 > /V3. With the origin of the position vector 
f at the center of the octahedron, the displacement vector 
u(r) of the corresponding pure shear is given by 

2 ^ ^ 

—u=(fa 2 )a 1 + (ra 1 )a2, (34) 
e c 

where e c is the shear angle of the gliding triangle mode 
of Fig. 6. 

From eq. (|34[) . one can calculate the shear. In Kittel's 
notation 30 (diagonal elements e xx — du x /dx, nondiago- 
nal elements e xy = du x /dy + duy/dx), one finds for the 
gliding triangle mode of Fig. 6 in units of e c v2/6 the 




FIG. 6: The octahedron and one of its twelve gliding triangle 
modes in cubic notation. Triangle ABC glides in the < 211 > 
direction, triangle DEF in the opposite direction. 
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diagonal elements e xx = —4, e yy — 2 and e zz = 2, the 
nondiagonal elements e xy = — 1, e xz = — 1 and e yz = 2. 
The elastic energy of the octahedron is 

E el = (C n - Ct2 + C AA )ve 2 c = (2C" + C iA )ve 2 c , (35) 

where C = (Cn — Ci2)/2 is one of the two shear moduli 
of the cubic crystal and C44 is the other. Equating both 
of them to the isotropic shear modulus G of the glass, 
we are back to the curvature of eq. © at the octahedral 
site. 

Two thirds of the elastic energy of eq. ([33)1 come from 
the diagonal elements, a dilatation in the y and in the 
z-direction accompanied by a factor of two larger con- 
traction in the x-direction. As it turns out, this special 
shear is a common feature of four of the twelve saddle 
point distortions, each coming from a different triangle 
plane. They only differ in the nondiagonal elements, the 
one from the (lll)-plane having e xy — 1, e xz = — 1 and 
e yz — —2, the one from the (lll)-plane having e xy = — 1, 
e xz = 1 and e yz = —2 and finally the one from the (111)- 
plane having e xy = 1, e xz = 1 and e yz = 2. This last one 
is closest to the first one, because the shear occurs in the 
same shear plane, with only a slight rotation of the shear 
directions. In the five-dimensional shear space, the lines 
connecting those two saddle points with the octahedron 
have an angle of only 38.9 degrees (Fig. 7). 

The closeness has consequences for the probability den- 
sity p s of section II. This was calculated integrating the 
partition function Z only over the connection line be- 
tween octahedron and saddle point, and then multiply- 
ing with the number of saddle points. In principle, one 
needs the full integration over the five-dimensional shear 



space. As long as the saddle points are well separated in 
phase space, the integration over the irrelevant degrees of 
freedom can be omitted^!. But if two saddle points are 
as close together as in Fig. 7, the integration over their 
neighborhood (along the dashed lines in Fig. 7) leads to 
double-counting. In a crude approximation, we counted 




saddle 1 



FIG. 7: The two neighboring saddle points with a strong 
common shear e common and a weak opposite shear e xy + e xz 
with respect to the undistorted octahedron. 

the four saddle points with the common diagonal shear 
elements as only two in eq. (|18|) (remember that the 
equation is only needed to demonstrate that the density 
of soft gliding triangle modes is comparable to the one of 
the sound waves). 
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